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Abstract — We report on terahertz frequency-domain spec- 
troscopy (THz-FDS) experiments in which we measure charge 
carrier dynamics and excitations of thin-film superconducting 
systems at low temperatures in the THz spectral range. The 
characteristics of the set-up and the experimental procedures are 
described comprehensively. We discuss the single-particle density 
of states and a theory of electrodynamic absorption and optical 
conductivity of conventional superconductors. We present the 
experimental performance of the setup at low temperatures for 
a broad spectral range from 3 to 38 cm"^ (0.1 - 1.1 THz) by the 
example of ultra-thin films of weakly disordered superconductors 
niobium nitride (NbN) and tantalum nitride (TaN) with different 
values of critical temperatures Tc. Furthermore, we analyze 
and interpret our experimental data within the framework of 
conventional Bardeen-Cooper-Schrieffer (BCS) theory of super- 
conductivity. By and large, we find the properties of our NbN 
and TaN thin films to be well described by the theory. Our results 
on NbN resemble tendencies towards anomalous behavior of the 
ratio 2A(0)/kBTc as a function of Tc. 

Index Terms — Frequency-domain THz spectroscopy, supercon- 
ducting thin films, BCS theory, density of states of a supercon- 
ductor, optical conductivity of superconductors, superconductor- 
insulator transition, TaN, NbN. 



I. Introduction 

THIN films of superconductors are of both technological 
and academic interest. Properties of interacting electron 
systems are governed by the electronic density of states 
(DOS) which changes upon reducing spatial dimensions. Thus, 
dimensional reduction 3D^2D may affect optical, electronic, 
and thermodynamic properties of the system. In superconduc- 
tivity, studying quasi-2D systems turned out to be a partic- 
ularly fruitful field and led to remarkable findings such as 
insulators featuring a superconducting gap [1], pseudogap in 
conventional s-wave superconductors [2], or extremely strong- 
coupling superconductivity in Kondo lattices |3|. Besides 
pure academic interest, thin-film superconductors play a key 
role in many state-of-the-art applications such as SQUIDs, 
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thermal and electrical switches using Josephson junctions, or 
microwave resonators [4]. Ultra- thin films of niobium nitride 
(NbN) and tantalum nitride (TaN) studied in this work are 
suited for a novel approach to build single-photon detectors 
|5|-|10|. Apart from applications, NbN has recently gained 
attention as a model system for the superconductor-insulator 
transition (SIT): Tunneling studies lUTIl reveal a gapped DOS 
slightly above Tc in highly disordered films as it has been 
found for titanium nitride (TiN) films [2J. This peculiar state 
resembles the well-established pseudogap in high-Tc cuprate 
superconductors. In case of NbN, however, Tc is about one 
order of magnitude higher compared to TiN. Thus, measure- 
ments on NbN in the vicinity of the SIT are much less 
demanding than in the case of TiN O, |[T? I. Consequently, 
results obtained for NbN might serve as a key to understand 
the still puzzling pseudogap state in high-Tc cuprates. 
THz time-resolved spectroscopy f\A\, THz frequency-domain 
spectroscopy (THz-FDS) 1 15 |, and far-IR Fourier- spectroscopy 
1 16 1 measurements on comparably thick films of NbN have 
been performed previously in a comparable spectral range. 
With our approach, however, we push experimental limits 
towards even lower frequencies and temperatures which allow 
us to study films with a higher degree of disorder and thus 
a smaller energy gap [llj. THz spectroscopy in the range 
from 0.03 to 1.5 THz (1 to 50cm~^ in the units common for 
spectroscopists) has proved to be a powerful tool to study 
many questions among thin film superconductors [17], since 
the superconducting energy gap 2 A of many compounds 
happens to fall in the corresponding energy range of the order 
of a millielectronvolt (~ 0.25 THz). Importantly, no contacts 
or surface structuring are required which might affect the sam- 
ple's properties, or complicate the analysis because of contact 
effects. The setup described below is a quasi-optical Mach- 
Zehnder interferometer that enables us to measure amplitude 
and phaseshift of THz radiation passing through a thin film 
on a plane-parallel substrate. From these experimental data, 
any complex optical function such as complex conductivity 
(J = (Ji + ia2 or the complex permittivity e = 6i + ie2 as well 
as free-electron parameters e.g., scattering rate i/r or plasma 
frequency ujp, can be calculated. 

This review paper is structured as follows. In section HI] we 
give a brief introduction in the BCS theory of superconduc- 
tivity focusing on the single-particle density of states and 
the electromagnetic absorption and the optical conductivity. 
Furthermore, in section |III] we describe our experimental 
techniques and the analysis approaches employed. Section JVl 
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features experimental data on thin film superconductors NbN 
and TaN, which not only play an interesting role in current 
research but also nicely demonstrate the agreement between 
theory and experiment. 

II. Electromagnetic Absorption and Optical 
Conductivity of a Superconductor 

In 1956, Cooper proposed an instability of a metaKs elec- 
tronic structure towards the condensation of free electrons 
into electron-electron pairs in the presence of a (arbitrary 
small) positive electron-electron net interaction [fTSi and, by 
this, laid the foundation for BCS theory of superconductivity 
|[T9l . By the exchange of a virtual phonon at sufficiently 
low temperatures, the repulsive Coulomb force between two 
(itinerant) electrons may be overcompensated lowering the 
energy of each electron by A(T), and lead to a bound electron- 
electron state, a Cooper pair, with the total binding energy 
2A(T). The phonon-mediated interaction is most effective in 
case of electrons with opposite momenta ki = k and k2 = — k 
and antiparallel spins si = 1/2 and S2 = — V2. Consequently, 
a Cooper pair has no net momentum and zero spin, thus being 
a boson. Obeying Bose-Einstein statistics, all Cooper pairs are 
allowed to occupy the same state which leads to the formation 
of a macroscopic phase-coherent condensate. Since the binding 
energy typically is of the order of a few millielectronvolts, the 
bound state is easily destroyed by thermal energy, and thus 
superconductivity commonly is limited to very low tempera- 
tures below a critical temperature Tc (Note: In what follows 
we will focus on conventional superconductors originating 
from metallic normal states. For reviews on unconventional 
superconductivity see e.g., ll2Ql - ll22ll ). 

A. Single-particle density of states 

The superconducting state features a number of intrigu- 
ing properties like non-dissipative transport current, perfect 
diamagnetism, and a peculiar single-particle density of states 
depending on temperature. At low temperatures in the super- 
conducting state, the step-like DOS of a metal features an 
energy gap of the size 2 A around the Fermi energy Ep and 
diverges at Ep ± A, see figure [T] The particular shape of 
the DOS sensitively depends on temperature. In the normal 
state, T > Tc, the metallic DOS is smoothly decreasing 
in the vicinity of Ep, and the energy states are occupied 
according to Fermi-Dirac distribution around the Ep, see 
figure [Ha). Because of the finite temperature, some states 
above Ef are occupied with thermal electrons while some 
states below Ef are empty. Upon cooling below Tc, the 
superconducting state evolves and an energy gap 2A(T) starts 
to open symmetrically around Ef relocating the states in the 
vicinity of Ef, see figure [Hb). The states which happen to 
fall in the energy gap are shifted above Ef -\- A(T) and 
below Ef — A(T) increasing the DOS at the edges of the 
energy gap. At Ef ± A(T), the DOS diverges and becomes 
infinite. The electrons which occupied the vanished states are 
bound into Cooper pairs and disappear from the single particle 
DOS when condensing to a mutual ground state. (In some 
models suited for e.g., tunneling phenomena. Cooper pairs are 
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Fig. 1. (Color online) Single-particle density states around the Fermi energy 
Ep = of 2i superconductor, (a) At T > Tc, the energy states are occupied 
according to Fermi-Dirac distribution leading to thermal population of some 
states above and some empty states below Ep. (b) At < T < Tc, a 
temperature-dependent energy gap 2A(T) opens around Ep and the DOS 
diverges at the edges of the gap. Due to finite temperatures, some states 
above the gap are occupied by thermally activated unpaired electrons. By 
absorption of photons of arbitrary energy, these thermal electrons contribute 
to (71. These transitions are marked by dashed arrows and correspond to 
integral A in equation [T] Electrons below the gap exclusively contribute for 
photon energies greater than 2A(T). These transitions are marked by solid 
arrows and correspond to integral B in equation [T] (c) At T = 0, no thermal 
electrons are apparent and 2A(T) reaches it's maximum. 



considered to be at But we employ a semiconductor-like 
picture, where Cooper pairs are not apparent). The majority 
of unpaired electrons remains below the gap, not contributing 
to the electrodynamic properties. A small number, however, 
is thermally activated to states above Ef -\- A(T). These 
thermally activated unpaired electrons (in the following just 
thermal electrons) play a key role for electrodynamics at finite 
temperatures. As the temperature approaches zero, 2A(T) 
rises to it's maximum and the number of thermal electrons 
decreases to zero, see figure [TJc). (In real systems, the DOS 
does not actually diverge but is smeared out due to the finite 
lifetime of the thermal electrons [|23l .) 
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B. Electromagnetic properties and optical conductivity 

The gapped and diverging DOS gives rise to significant 
changes in the electromagnetic absorption and the optical 
conductivity for frequencies of the same order as the gap 
frequency uOg = 2A/h, with h = h/27r the reduced Planck 
constant. Considering typical values of the energy gap, these 
changes are likely to be found in the far-IR, THz and mi- 
crowave spectral ranges. The optical conductivity of a super- 
conductor was first derived by Mattis and Bardeen [24] and 
holds for superconductors, where the mean free path i oc r, 
with r the relaxation time, is by far superior to other typical 
length scales. The more general formalism for arbitrary values 
of T is first presented in [25] , and casts into the equation (with 
Ef=0) 
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1 1 = tanh 



range [A(T),oo) since thermal electrons occupy states up to 
arbitrary high energies. Integral A covers the energy range 
[A(T), huj-\-A{T)], however, concerning ai it yields non-zero 
values only up to huj — A(T).) For frequencies smaller than 
ujg, however, the photon energy becomes insufficient to excite 
electrons across the gap, and consequently the electromagnetic 
response is completely governed by intraband transitions of 
thermal electrons. The number of possible transitions by 
photon absorption increases with decreasing frequency due 
to the diverging DOS. Consequently, ai rises strongly as 
the frequency approaches zero. At a certain frequency, ai 
exceeds the normal state conductivity adc- This becomes 
obvious from spectral weight arguments, since the area 
beneath ai{uj) is constant for all temperatures according to 
Ferrell-Tinkham-Glover sum rule [27], [28 1. The missing 
area beneath ai{uj) in the superconducting state is relocated 
to both the (5-function at zero frequency and to the low- 
frequency range, where ai rises above the normal state 
conductivity. For uj = 0, the contribution of thermal 
electrons merges with the infinitely high S-likc conductivity 
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and ks the Boltzmann constant. (For numeric calculations, a 
more convenient version of equation [T] is given in [25 1.) Real 
and imaginary parts of the complex conductivity follow from 
equation [T] with ai = Re((j) and (72 = Im((j). 
The first integral hdE in equation [T] (in the 

following labeled integral A) accounts for the electromagnetic 
response of the Cooper pairs, while the second one 
/a(t) ^2^^ (integral B) accounts for the thermal electrons. 
Typical graphs for cri{uj) and (72 (cj) calculated using equation 
[T] are shown in figure |2] Besides the total response taking 
into account both integrals A and P, figure [2] disentangles 
the contributions resulting from Cooper pairs and thermal 
electrons. For finite temperatures below Tc and frequencies 
higher than ujg, electrons are activated across the energy 
gap by absorbing photons, see figure [T]3), thus giving rise 
to a finite (7i (This process corresponds to the breakup 
of a Cooper pair, adding two quasip articles - a quantum 
superposition of electron and hole wavefunctions e.g., [26] - 
to the single-particle DOS). With increasing photon energy, 
the range from where electrons may be excited across the 
gap becomes greater enabling enhanced absorption and 
consequently rising cri. The contribution of thermal electrons 
is finite, however, inferior to the previously discussed 
contribution. This is mainly due to the comparably small 
amount of thermal electrons and the small DOS at energies 
higher than A(T) + fiuog. (Note: integral B covers the energy 



of the superconducting condensate at zero frequency. (Note: 
equation [T] does not include the ^-function contribution to 
ai.) 

While the real part cfi of the optical conductivity features 
both considerable contributions from integrals A and P, the 
finite imaginary part (72 is mainly due to the superconducting 
condensate (at sufficiently low temperatures and frequencies), 
whose conductivity is derived via the 1. London equation 

\m 

djs{uj,t) nsC^ 



dt 



m 



-E{uj,t) 



m 



E(uj,t')dt' 



(2) 



(with Us the superfluid density, e the elementary charge, and m 
the electron mass) relating the time derivative of the (local) su- 
percurrent density js(^, to an external perturbation E{uo^t). 
In case of an electromagnetic wave, E{uj^t) = Eq exp(— icjt), 
equation |2] then becomes 

o 

1 77/ 6 

js{cj,t) = —E{^,t) (3) 



UJ m 



Equation [3] has the same structure as Ohm's law, and thus we 
interpret the proportionality factor ^ between js (cj, t) and 
E{uj,t) as the conductivity of the superconducting condensate. 
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Frequency (2A(T)/M) 



Fig. 2. (Color online)(a) Real and (b) imaginary parts of the complex optical 
conductivity (solid line) of a superconductor as a function of frequency based 
on equation [T] for T = 0.8Tc. Dashed lines disentangle the contributions 
of both the thermal electrons (Integral B in equation [T) and photon-activated 
electrons (Integral A in equation [TJ to the conductivity. Insets of (a) and (b) 
show (71 and (72 for different temperatures (T/Tc = 0.9, 0.8, 0.7, 0.6, 0.5, 0.4, 
0.3, 0.2, and 0.1) and scattering rates (i/r = 40, 20, 4, 2, 0.8, 0.4, and 0.04 
x2A(0)), respectively. (Used parameters 2A(0) = 25 cm^^, Tc = 10 K, 
(Jdc = 10661^-icm-i) 



which is purely imaginary and thus only contributes to (T2. 
Upon increasing frequency, (72 decays as i/u;. Consequently, 
the strongest contribution is at energies well below 2A(T), 
much lower than addressed in our study, i.e. at microwave 
frequencies ll30ll , ll3T1l . Kramers-Kronig analysis underlines the 
above qualitative picture yielding the same result [29|. 
Integral A yields finite values up to huj + A(T) which calls 
for the integration range given in equation [T] For frequencies 
larger than ujg, Cooper pairs are destroyed, and thus the overall 
out-of-phase response changes. This is signaled by a weak 
kink of (72 at the energy gap, see figure |2] In analogy to an 
electric LC circuit, the frequency dependence of (72 is often 
interpreted in terms of frequency dependent kinetic inductance 
I32II . I33I . The impedance of the circuit rises towards higher 
frequencies leading to a reduced conductivity. A finite scat- 
tering rate i/r accounted for in equation \T\ further influences 
the behavior of the complex conductivity. Upon increasing 1 /r 
(going from clean-limit to dirty-limit superconductors), the rise 
in (7i for uj > ujg is gradually suppressed, and (72 features a 
strongly pronounced kink at ujg for scattering rates comparable 
to the energy gap, see insets of figure Oa) and (b). 



III. Experiment and data processing 
A. THz frequency -domain spectroscopy 

Our THz studies are based on the frequency-domain spec- 
troscopy technique. Coherent, monochromatic and continuous 
THz radiation is generated by frequency-tunable backward- 
wave-oscillators (BWO) that provide frequencies from 1 to 
47cm~^ (341. This spectral range is covered by a number 
of different BWO radiation sources (in the following just 
sources), see bottom of figure [3] The power of the radi- 
ation strongly depends on the frequency and may be as 
large as lOOmW for low-frequency and significantly smaller 
(0.1 mW) for high-frequency sources. As a consequence, the 
low-frequency radiation power might be too strong and re- 
quires use of attenuators to avoid sample heating or detector 
overload. The great benefit of FDS compared to time-domain 
spectroscopy (TDS) is to process frequencies step by step with 
a stability and resolution exceeding Au/uj = 10~^ depending 
on the BWO power supply. 

While TDS excites all available energies simultaneously and 
has to struggle with minor signal strength, FDS can directly 
probe narrow-band excitations and resolve detailed line shapes. 
The high output power leads to high signal-to-noise ratio that 
can reach values up to 10^. Frequency doublers and triplers 
(Virginia Diodes) that allow extending the frequency range of a 
separate source can effectively be used. Loss of the radiation 
intensity is easily compensated by use of a ^He bolometer 
(see below). After the radiation is generated by the BWO, it 
is shaped by a Teflon, quartz, or polyethylene lens and guided 
by aluminum mirrors and wire grids, see top of figure |3] which 
act as beam splitters depending on polarization. Other lenses 
are used to focus the radiation on the sample under study and 
the detector. Attenuators, polarizer, chopper, and analyzer are 
tilted with respect to the beam propagation axis to suppress 
standing waves between these parts. The quasi-optical part of 
the experiment is a Mach-Zehnder interferometer which allows 
us to measure both amplitude and phaseshift of THz radiation 
passing through a sample. The radiation is detected by either a 
Golay cell or a "^He-cooled bolometer (Infrared Laboratories), 
depending on signal strength. For transmission measurements, 
the radiation is mechanically chopped to make it detectable 
using a lock-in technique. 

Determining the spectra of optical parameters of the plane- 
parallel sample involves measurements of the transmission 
and phaseshift spectra from which the spectra of real and 
imaginary parts of complex permittivity, conductivity, etc., are 
calculated directly and without use of any additional analysis 
(such as Kramers-Kronig analysis). 

Recording a transmission spectrum is performed with the 
second arm of the interferometer blocked. It consists of two 
separate steps. First, the signal intensity versus frequency is 
recorded with no sample in the beam path. This transmis- 
sion spectrum is regarded 100% transmission and used as 
measurement calibration. Second, the transmission spectrum 
is recorded with the sample in the beam path and the absolute 
transmission spectrum is obtained as a result of division of the 
two correspondent data arrays. In this way we can account 
for the frequency-dependent output power of the radiation 
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Fig. 3. (Color online) Top: Simplified sketch of the utilized Mach-Zehnder 
interferometer. Attenuators, polarizer, chopper, and analyzer are tilted with 
respect to the beam propagation axis to suppress standing waves between 
these parts. Measuring the phaseshift requires both arms of the interferometer, 
whereas the second arm is blocked for transmission measurements. Bottom: 
(a) BWO output power of several sources versus frequency. The range from 
1 to 47 cm~ ^ is covered entirely, (b) The available spectral range of a source 
can be increased using passive frequency doubler and tripler (dashed lines). 
The signal intensity in the doubled/tripled spectral range, however, is 3 to 4 
orders of magnitude smaller and may require use of a ^He-cooled bolometer 
for detection. 



passing through the measurement channel. Typical data (300 
K) of a 0.78 mm thick substrate of weakly doped silicon 
(n = 3.5 X lO^^cm"^) are shown in figure IHa) [12J. The 
oscillation pattern is due to multiple reflections inside the 
sample which acts as a Fabry-Perot resonator. 
To measure phaseshift spectra, the second arm has to be 
unblocked. Again, the measurement consists of two separate 
steps: calibration without and measurement with the sample 
in the beam path. In order to obtain reliable data, one has to 
guarantee optimal interference of the two beams at the ana- 
lyzer. This is provided by aligning the spectrometer elements 
and achieving maximal ratio between minimal (destructive 
interference) and maximal (constructive interference) signal 
strength. The interference signal is detected using a lock- 
in technique. However, the signal is modulated not using 
the chopper (amplitude modulation) but the oscillating mirror 
(phase modulation), see figure O During a measurement 
process, the movable mirror is automatically kept in a position 
corresponding to zero-order destructive interference (minimal 
signal) while the frequency range is swept. The mirror po- 
sition versus frequency is recorded when the sample is out 
of (calibration scan) and in (measurement) the optical path, 
and the phaseshift spectrum is calculated from the difference 
between the two spectra. Typical phaseshift data taken on the 
same 0.78 mm thick substrate of doped silicon is provided 
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Fig. 4. (Color online) Typical raw data of (a) transmission and (b) 
relative phaseshift versus frequency together with simultaneously performed 
fits according to equations |4] and |5] taking into account Drude behavior given 
by equation [21 The sample is a ~ 0.78 mm thick substrate of weakly doped 
silicon at room temperature. Note: the phaseshift is the measured quantity, 
however, it is common to work with phaseshift divided by frequency. This 
will be done throughout the entire work. 



in figure Ob). (Note: the measured quantity is phaseshift (j) 
(rad), however, in this paper we will always refer to relative 
phaseshift (rad/cm~^) with u = uj/{2ttc) and c the speed 
of light in vacuum). 

Optical anisotropy is an issue that has to be taken into 
account carefully. E.g., the Fabry-Perot resonances can become 
more complex |35|. Many common substrates for thin-film 
deposition are birefringent such as AI2O3 or NdGaOs. To align 
the polarization direction E of the radiation to sample's axes 
of interest, one can adjust either the sample or the polarization 
angle of the 1. beam splitter (followed by corresponding 
turning of the 2. beam splitter, setting it to 90° relative to 
the 1. beam splitter). In case of optical anisotropy of both 
substrate and film and misaligned optical axes, the optical set- 
up and analysis has to be customized 



B. Data analysis of single-layer systems 

Transmission and phaseshift data are fitted by well-known 
Fresnel equations for multiple reflections at media boundaries. 
For a single layer system, the transmission T and phaseshift 
6 read 129,1 



T 



[{I - Rf + ARsm^ if^ exp{-a(i} 
(1 - i?exp{-ad)}2 + 4i?exp{-ad} sin^ (/3 + ip) 



(4) 



2iind 



arctan 



I (2 + n) n 



arctan 



Rexp{— ad} sin^ -\- (p) 
1 — R ex.p{—ad} cos^ {13 -\- (p) ^ 



(5) 
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where = (l-2n+n^+/c^)(2n+n^+/c^)~^ is the reflectivity, 
= arctanj— 2/c(l— n^ — is the phase change upon re- 
flection at an interface between the two media, /3 = 27rnd/Xo 
is the angle by which the phase of the radiation is changed 
upon TraveUing through a medium of thickness d for a given 
wavelength Aq and a = Ank/Xo is the power absorption 
coefficient. The quantities n, /c are real and imaginary parts 
of the complex refractive index >/I = n = n + i/c in the given 
case of a non-magnetic material. Equations |4] and [5] do not 
account for any frequency dependence for the optical constants 
n and k, which typically depend on frequency and reveal the 
properties of interest in this study. To include resonant absorp- 
tion (e.g. Lorentz oscillators) or conductivity (Drude behavior), 
corresponding terms are used. For example, to describe the 
free carrier response within the Drude conductivity model we 
can introduce complex permittivity e = ei + ie2, and e.g., in 
case of a simple metal the Drude formula [.29J 

, ,2 



Frequency (THz) 
0.16 0.20 



e2(^) 



1 



(6) 



(7) 



ur -\- r 

or rewritten in terms of the complex optical conductivity 

(7H=i^(eoo-e(a;)) 

cri(cj) 



CTdc 



1 



a2{uj) = (Jdc 



1 



(8) 



(9) 



where ujp is the plasma frequency, cfdc = ^^V^^ the dc 
conductivity, eoo is the value of the dielectric constant for 
frequencies much higher than addressed in this study, and 
i/r is the relaxation rate. The solid lines in figure |4] are 
fits based upon equations H] and [5] The fitting parameters 
n and k are expressed in terms of ei,e2,cri and a2 using 
well-known expressions, [29 1, and, by this, the frequency 
dependence according to Drude theory enters the transmission 
and phaseshift fits via equations [6l|9l Clearly, this yields 
an excellent description of the experimental data, i.e. this 
silicon substrate is a text book representative of optical Drude 
behavior i37ll-ll39ll. 

Another analysis approach is to fit the Fabry-Perot transmis- 
sion peaks separately, with two advantages: first, the noise 
is least pronounced in the vicinity of the Fabry-Perot max- 
ima, and the largest sensitivity is reached here. Second, this 
routine provides more freedom in the case of materials with 
unknown properties where no explicit function for frequency 
dependence is known. To do so, each Fabry-Perot maximum 
at some given frequency cjq is fitted separately and yields a 
set of material parameters at uq. Processing each Fabry-Perot 
transmission peak in this way eventually yields the desired 
frequency-dependent optical functions. 

C. Data analysis of multi-layer systems 

Equations (H and \5\ hold for a single layer. For a two (or 
more)- layer system, however, another boundary has to be 
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Cryostat 2 (tilted windows) ' 
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Fig. 5. (Color online) Raw transmission (a) and relative phaseshift (b) versus 
frequency measured with parallel-outer- window cryostat "1" (filled circles) 
and tilted-outer-window cryostat "2" (open circles) at two temperatures 4.2 
and 12 K (sample: NbN, Tc = 11.8 K). Towards low frequencies, standing 
waves are strongly pronounced for cryostat "1". With increasing frequency the 
standing waves become less pronounced. Simplified sketches of the window 
arrangement in top- and side-view positions of both cryostats are shown in 
(b). 



accounted for. To disentangle properties of the two media (e.g. 
film on a substrate) the two experimental quantities (trans- 
mission and phaseshift) have to be measured. The response 
functions of the film are calculated from the two correspondent 
spectra. This requires detailed knowledge of the substrate 
parameters, which have to be determined beforehand in a 
separate experiment. The Fresnel equations are extended by 
terms describing another boundary reflection and are found 
elsewhere e.g., ll29ll or BOl . To improve the interaction of 
probing radiation with the material under study, a resonant 
Fabry-Perot technique has been developed where two speci- 
mens of the thin film under study are mounted face-to-face 
realizing a resonant cavity BTl . This technique is useful in 
case of superconducting films, where it can be difficult to 
obtain the dynamical conductivity ai with sufficient accuracy 
due to huge positive values of (72. Fabrication of such a 
resonant cavity, however, remains a challenging task [41]. 

D. Optical cryostats for THz spectroscopy and standing-wave 
problem 

To access low temperatures we employ two different home- 
built ^He-bath cryostats labeled "1" and "2" here. The most 
notable difference among these cryostats is the arrangement 
of the optical inner and outer windows. Diffraction and stand- 
ing waves distorting the measured spectra become a severe 
problem especially at low frequencies below 10cm~^, since 
the THz radiation wavelength becomes of the same order as 
typical setup dimensions, e.g. the distances between inner 
windows of the cryostat, polarizers, lenses, etc.. To reduce 
this problem, the windows are made of thin Mylar (~ 20 fim) 
foils which are slightly curved due to inner vacuum. The outer 
windows of cryostat "1" are oriented perpendicularly to the 
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Fig. 6. (Color online) Sheet dc resistance versus temperature for different 
films of (a) TaN and (b) NbN. Different values of Tc are obtained for different 
deposition currents while sputtering. We estimate the loffe-Regel parameter 
kpi to vary between 4.5 and 5.5 for the NbN films under study far off the 
superconductor-insulator transition. 

beam still leading to interference which results in multiple 
reflections and noticeable standing waves. Cryostat "2" was 
designed with both tilted inner and outer windows (tilting 
angle ~10°), see inset of figureOb). Transmission and relative 
phaseshift spectra taken with both cryostat "1" and "2" are 
displayed in figure [5] for comparison. Clearly, standing waves 
are much less pronounced on cryostat "2". Obviously, the tilted 
windows have to be larger than the non-tilted ones to not 
reduce the beam cross section. Consequently, the amount of 
thermal radiation entering the sample chamber is bigger which 
introduces additional heating, complicating thermal control 
at temperatures below the boiling point of ^He at ambient 
pressure. 

Another common approach to reduce the standing wave pattern 
in the measured spectra is modulating the electron- acceleration 
voltage (0.5 to 6kV) applied to the BWO, i.e. modulating the 
frequency of the THz radiation making it less monochromatic. 
Narrow-band features of the spectra might be washed out 
by this frequency averaging, and therefore the high- voltage- 
modulation amplitude Vp-p (0 to 50 V) has to be adjusted 
to find an optimum between reduced standing waves and 
frequency fidelity. 

IV. Experiment 

The thin films under study were deposited on 10 mm x 
10 mm R-plane-cut birefringent sapphire substrates using re- 
active dc-magnetron sputtering. The NbN and TaN films have 
a thickness of 3.9 to 4.7 nm and 4.9 to 5.1 nm, respectively. 
Depending on the deposition current Id while sputtering, the 
critical temperature is lowered in case of NbN from 12.2 K 
(Id = 100 mA) to 8.5 K (200mA) and in case of TaN from 
9.7 K (70mA) to 8.6 K (120mA), see figured Comparison 
with NbN film deposited in the identical way [42] and [llj, 
allows us to estimate the amount of disorder via the loffe- 
Regel parameter kp^ (Fermi wave vector x electron mean free 
path) that varies between 4.5 and 5.5 for our NbN films. Thus, 
we consider our samples to be only weakly disordered and far 
off the SIT. For more details on the samples and deposition 
process see BSl and |6|. We investigate several NbN and 
TaN samples at temperatures above and below Tc in the 
frequency range 3.4 to 38 cm~^. To cover this frequency range, 
5 different sources were used. Starting at low frequencies, the 



output power of the first 4 sources was high enough to be 
detected by a Golay cell. In contrast, to achieve a reasonable 
signal-to-noise ratio with the high-frequency source, a ^He- 
cooled bolometer was necessary. Transmission and phaseshift 
measurements were carried out on both main axes of the 
birefringent sapphire substrate. The samples were mounted 
into the optical cryostats in such a way that radiation hits 
the substrate first and subsequently the film. 

V. Results 

Comprehensive low temperature THz measurements were 
performed on all TaN samples using cryostat "2" exclusively. 
In case of the NbN samples, only measurements at 4.2 K were 
carried out with "2" otherwise "1". Due to similar results we 
discuss only one TaN and one NbN sample on behalf of the 
others. Raw data of transmission and relative phaseshift versus 
frequency above and below Tc are shown in figure |7] 



Frequency (THz) 
0.2 0.4 0.6 0.8 1.0 




10 20 30 



Frequency (cm"^) 

Fig. 7. (Color online) Raw transmission of (a) NbN and (c) TaN and relative 
phaseshift of (b) NbN and (d) TaN versus frequency at two temperatures, 4.2 K 
and in the normal state. Solid lines are transmission and relative phaseshift fits 
with the optical conductivity given by equations |6] and [7] in the normal state, 
and equation [T] at 4.2 K. The vertical dashed lines mark the gap frequency 
ug = ujg/(27Tc) ^ 28cm~i (NbN) and Ug = ujg/(27Tc) ^ 21.7 cm~i 
(TaN). 
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In general, raw data of both NbN and TaN are very 
similar. In the normal state, the overall Fabry-Perot pattern 
of the transmission coefficient and phaseshift does not show a 
frequency dependence. This indicates a constant dispersionless 
complex conductivity of the film, similar to those of a simple 
metal for frequencies much lower than the scattering rate. The 
transmission reaches 45 to 50% signaling a comparably poor 
metallic conductivity. When the samples are cooled below 
Tc, transmission and relative phaseshift spectra change signifi- 
cantly. The amplitudes of the Fabry-Perot transmission peaks 
exceed the normal-state values for all studied frequencies. 
Starting at low frequencies, the maxima in transmission in 
NbN increase from 0.7 to 1 at about 24cm~^ and slowly de- 
crease to 0.9 at high frequencies. Similar behavior is apparent 
for TaN, however, the maxima in transmission decline more 
rapidly towards high frequencies. Towards low frequencies, 
relative phaseshift of both NbN and TaN is strongly reduced 
compared to the normal state. Amplitude and periodicity of 
T and ^/cj vary with temperature and frequency signaling 
a considerable frequency dependence of the complex optical 
conductivity. 
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Fig. 8. (Color online) (a) Energy gap 2 A versus temperature of a NbN 
sample for both main axes of the substrate. Solid lines are fits according to 
BCS theory, equation [To] Results on both main axes yield very similar values 
of 2A(T). The error bar shown is representative for all data, (b) 2 A at 4.2 K 
and (c) ratios 2A(0) /A^^Tc versus Tc of all studied NbN films. 



VI. Analysis and Discussion of NbN 

All raw transmission and relative phaseshift spectra for 
the superconducting state measurments are fitted taking into 
account the frequency dependence of the conductivity given 
by equation [T] For each temperature, the fits for transmission 
and relative phaseshift are performed simultaneously with 
2 A being the only fit parameter. All free-electron parameters 
in equation [T] (scattering rate i/r and dc conductivity a^c) 
are obtained from the normal- state measurements directly 
above Tc, which are fitted to Drude theory and taken as 
constant. The critical temperature Tc is obtained from transport 
measurements, see figure |6] BCS fits for both transmission 
and relative phaseshift at 4.2 K are shown in figure [71 Clearly, 
theory and experiment match very well for frequencies both 
above and below the energy gap. The same holds for the other 
temperatures, which allows to extract 2 A for each temperature 
precisely. Figure [8] displays the temperature dependent energy 
gap 2A(T) of the sample with lowest disorder and highest 
Tc. Filled and open symbols refer to different main axes of 
the substrate and coincide within the range of error (52a = 
0.9 meV. Data for each main axis is fitted separately with the 
BCS expression |l441 



A(T) 



tanh 



Tc A(T) 1 



(10) 



A(0) i T A(0) J 

For 0.3Tc < T < 0.8Tc the measured energy gap is slightly 
above the BCS curve and for 0.8Tc < T < Tc below. The 
low-temperature part displays a similar behavior as it has 
been observed in microwave experiments [45], where 2A(T) 
does not decrease as fast as in BCS theory as T Tc. 
In our measurement, however, this trend is not true for 
temperatures close to Tc and our findings do not signal a 
gapped normal state. Upon decreasing Tc, the energy gap 
(at 4.2 K) is reduced from about 4.5 to 3.3 me V, see inset 
of figure m Consequently, 2A(T) depends on the amount of 
disorder and is greatest for low-disorder films. Furthermore, 



we do not see optical signatures of superconductivity in the 
normal state among the higher disordered samples. At 4.2 K, 
2A(T) is reduced as Tc decreases, which underlines previous 
findings from tunneling studies ifTTIl . see inset of figure [8] 
We calculate the ratios 2A(0) //c^Tc from our data assuming 
2A(4. 2 iiT) ^ 2A(0), and find a similar anomalous behavior 
as in lITTIl : upon increasing disorder (decreasing Tc), the ratio 
2A(0) /ksTc is first reduced from about 4.1, indicating strong- 
coupling superconductivity, exhibits a minimum and becomes 
larger again (about 4.5) towards the highest-disorder sample 
studied. Our overall findings, however, favor conventional BCS 
superconductivity and do not feature any superconducting 
properties that require interpretation in terms of the SIT. This 
is expected to change as the film thickness is reduced. 

VII. Analysis and Discussion on TaN 

The experimental data for TaN is processed in the same 
way as described for NbN. We find our experimental data in 
perfect agreement with weak-coupling BCS theory, see figures 
[7] and [9] and [TOl representative for all samples under study. 
From the BCS fits for all temperatures below Tc we obtain 
the temperature-dependent energy gap, which is displayed 
in figure |9l Considering responses measured along the main 
axes individually, theory and experiment match nicely. Within 
our range of error 82 a = 0.9 me V, results on both axes 
(denoted by "+" and "-" ) are the same. We attribute the 
fairly constant difference for both axes to a small systematic 
experimental error aligning the substrate's axes and the THz 
polarization. The inset of figure [9] displays 2 A at 4.2 K for 
all studied samples versus Tc. While 2 A tends to decrease as 
Tc is reduced in case of NbN, we do not see a significant 
dependence for TaN. Since the ratio between highest and 
lowest critical temperatures t^^^ /T^^^ is 1.13 and 1.44 for 
TaN and NbN, respectively, we can only access a smaller range 
here which makes it difficult to establish similar conclusions 
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Fig. 9. (Color online) (a) Energy gap 2 A versus temperature for both main 
axes of the substrate (denoted by "+" and "-" ). Solid lines are fits according 
to BCS theory, equation [Tol Considering each axis separately, experiment and 
theory match nicely. The error bar shown is representative for all data, (b) 
2 A at 4.2 K and (c) ratios 2A(0) //c^Tc versus Tc of all studied TaN films 



from our TaN experiments. We calculate the ratio 2A(0)//cbTc 
from our experimental data to be closely scattered around 
the weak-coupling prediction 3.52 lfT9ll , which underlines pure 
weak-coupling BCS superconductivity in this range of disorder 
and critical temperatures. 

We also analyze the raw data with single-peak fits to ex- 
tract the frequency dependence of the complex conductivity 
d{uj) = cri(cj) + \(T2{(jo) independently from a predefined 
model. For this, a narrow frequency range of 0.5 to 1 cm~^ is 
chosen for fitting each single Fabry-Perot peak separately with 
ei and 62 as fitting parameters. Subsequently, cri and cf2 are 
calculated via e = eoo + i^ircr/uj for each peak frequency with 
600 = 1. Figure [TOl displays the conductivity for 4.2 and 7 K 
well below Tc = 9.7 K and at 20 K in the normal state. Closed 
and open symbols are for the different main axes. Data points 
corresponding to different main axes are shifted in frequency. 
The reason is that ei of the birefringent substrate is slightly 
different for both axes leading to a shift in peak position. The 
solid lines are BCS predictions based on the raw-data fits. 
Within the range of error, theory and experiment (for both 
main axes) match reasonably. 2A(T) is signaled by the kink in 
(Ji and is clearly resolved within the experimental frequency 
range. At 20 K, cfi does not depend on the frequency and 
corresponds to the constant value ai = cfdc = 4580 ft~^cm~^ 
extracted from Drude fits. The studied frequency range is well 
below the scattering rate i/r and, thus, (72 is immeasurably 
small in the normal state at 20 K. In the superconducting state, 
however, (72 rises notably and approaches the zero-temperature 
prediction (72 (T = 0) ex uj'^ for huj < 2A(T) , which 
is consistent with the ^-function-like conductivity ai at zero 
frequency. This trend in (72 is also seen in the superconducting 
penetration depth A l • The inset of figure \W[h) displays A l 
versus frequency for both main axes at 4.2 and 7 K calculated 



(11) 



For frequencies Tiuj < 2A(T) and a (^-function-like (7i, 
Al does not significantly depend on frequency. Despite the 
finite temperatures, such frequency-independent behavior is 
supported by the data in figure [TOl At 4.2 K, Al remains fairly 
constant for frequencies lower than 24cm~^ and approaches 
~ 500 nm for a; ^ 0. At 7 K, similar behavior is apparent for 
frequencies lower than 15cm~^, and Al approaches ~ 600 
nm. The zero-frequency values of Al are very similar to ones 
obtained from previous transport measurements [46 1. By and 
large, we find the charge carrier dynamics of our TaN films 
to be perfectly described within weak-coupling BCS theory. 
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Fig. 10. (Color online) Real part (a) and imaginary part (b) of the complex 
conductivity versus frequency for 4.2 and 7K well below Tc and 20 K in 
the normal state for both main axes of a TaN sample. Solid and dashed 
lines are predictions calculated directly from raw data fits for both axes 
individually. The shift in frequency is caused by different peak positions 
due to birefringence. The kink in ai signaling the energy gap is resolved 
clearly. The error bar shown is representative for all data. The inset shows the 
penetration depth as a function of frequency. A/, is frequency independent for 
TiLo < 2A(T). 



VIII. Summary 

We have described the performance and characteristics of 
our frequency-domain THz spectroscopic experimental tech- 
niques. Furthermore, we give a comprehensive overview of 
the analyses we employ and the theoretical foundation on 
which they are based. We studied charge carrier dynamics 
and superconducting properties of NbN and TaN thin films 
with different values of Tc at low temperatures above and 
below Tc in the THz frequency range from 0.09 to 1.1 THz 
(3 to 38cm~^). Coherent and continuous THz radiation was 
generated by frequency-tunable high power backward wave 
oscillators covering the frequency range below and above 
the energy gap 2A(T). We fit our experimental data of 
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transmission and phaseshift to a combination of Fresnel equa- 
tions, BCS- and Drude theory and describe superconducting 
properties and charge carrier dynamics. The results for both 
NbN and TaN are very well described within the framework of 
BCS theory. The agreement between our results and previous 
measurements on NbN from another group confirm that our 
experimental approach is perfectly suited to study correlated 
electron phenomena in thin film systems. The dependence of 
the ratio 2A{0)/kBTc and 2A on Tc we obtain for our NbN 
films resemble anomalous behavior as it has been reported 
previously, whereas the ratios for our TaN films do not 
significantly deviate from the weak-coupling BCS prediction. 
Absolute values of 2A and the ratio 2A(0)/A:bTc for NbN 
are in good agreement with previous publications, and thus 
underline the reliability and accuracy of our experiment and 
analyses. This also endorses our results on TaN, a compound 
that so far has not been studied in detail concerning THz 
charge carrier dynamics. 
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